arXiv:1507.04358vl [astro-ph.GA] 15Jul2015 


Draft version July 17, 2015 

Preprint typeset using DT^X style emulateapj v. 5/2/11 


MEASURING THE MASS OF THE CENTRAL BLACK HOLE IN THE BULGELESS GALAXY NGC 4395 

FROM GAS DYNAMICAL MODELING 

Mark den Brok 1 7 , Anil C. Seth 1 , Aaron J. Barth 2 , Daniel J. Carson 2 , Nadine Neumayer 3 , Michele Cappellari 4 , 
Victor P. Debattista 5 , Luis C. Ho 6,7 , Carol E. Hood 8 9 10 , and Richard M. McDermid 9,10 

Draft version July 17, 2015 

ABSTRACT 

NGC 4395 is a bulgeless spiral galaxy, harboring one of the nearest known type 1 Seyfert nuclei. 
Although there is no consensus on the mass of its central engine, several estimates suggest it to be one 
of the lightest massive black holes (MBHs) known. We present the first direct dynamical measurement 
of the mass of this MBH from a combination of two-dimensional gas kinematic data, obtained with the 
adaptive optics assisted near infrared integral field spectrograph Gemini/NIFS, and high-resolution 
multiband photometric data from Hubble Space Telescope's Wide Field Camera 3 {HST/ WFC3). We 
use the photometric data to model the shape and stellar mass-to-light ratio (M/L) of the nuclear star 
cluster. From the Gemini/NIFS observations, we derive the kinematics of warm molecular hydrogen 
gas as traced by emission through the H 2 1-0 S(l) transition. These kinematics show a clear rotational 
signal, with a position angle orthogonal to NGC 4395’s radio jet. Our best fitting tilted ring models of 
the kinematics of the molecular hydrogen gas contain a black hole with mass M = 4l® x 10 5 M0(3er 
uncertainties) embedded in a nuclear star cluster of mass M = 2 x 10 6 Mg. Our black hole mass 
measurement is in excellent agreement with the reverberation mapping mass estimate of Peterson et 
al. (2005), but shows some tension with other mass measurement methods based on accretion signals. 

Subject headings: galaxies: spiral - — galaxies: individual (NGC 4395) — galaxies: nuclei — galaxies: 
kinematics and dynamics — galaxies: structure 


1. INTRODUCTION 

It is suspected that almost all high-mass galaxies 
harbor a supermassive black hole at their center (e.g . 
IKormendv fc Richstonel 119951 : iKormendv fo Hoi 120131 1. 
There are two popular scenarios to form such massive 
black holes. In the first scenario, black holes are cre¬ 
ated as stellar remnan ts, most likely from hea v y Pop¬ 
ulation III stars ( e.g. iBromm fe Larsoiil 12004 Ig lover! 
2005; Br omml 1201311 . The second option is a scenario in 
which primordial black hole seeds form from direct col¬ 
lapse of_gas_clouds_Jeadingto 10 1 2 3 4 -10 5 black holes 
(e.g. lHaehnelt fe Reesl 119931 : iKoushiappas et ahl 12004 
iBegelman et al.ll2006t) . Both scenarios would require a 
period of rapid growth for the seed black hole in order 
to produce the massive quasars found at high redshifts 
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dFan et, al.l 12004 l.liang et al.l 120091 : iWillott et al.l 120091 : 
iVenemans et, al.ll2013fl. 

Due to subsequent growth of these seeds, both black 
hole formation scenarios predict high present day oc¬ 
cupation fraction and the same black hole mass func¬ 
tions in high-mass galaxies. However, the direct col¬ 
lapse scenario predicts a much lower occupation frac- 
tion in low-mass galaxies than th e remnant scenario 
(IVolonteri et al.! 12008b [G reen llpl. A census of low- 
mass supermassive black holes in low-mass galaxies can 
thus provide important constraints on the formation 
mechanisms of supermassive black holes in the early Uni¬ 
verse. For high-mass galaxies, the mass of the MBH 
appears to correlate tightly with the properties of the 
bulge. Specifically, there is a correl ation with the lu¬ 
minosity and mass of the bu lge (IKormendv fe Richstonel 
119951: [MagoiTia.li et al.l Il998h . but the strongest correla¬ 
tion was found to be between black hole mass and veloc- 
ity dispersion of the stars in the bulge, the M.-q r elation 
(|Ferrarese fc Merritt! 120001 : IGebhardt et al.ll2000H . The 
tightness of this relation implies some physical connec¬ 
tion between the evolution of bulges and MBHs. Com¬ 
bined with the lack of a massive black hole in M33, this 
led to the suggestio n that black holes may be absent from 
pure disk systems (|Merritt et al.1120011 : IGebhardt et al.1 
1200 If ). The presence of a black hole in NGC4395 defies 
thi s expectati on. The galaxy is classified as a Sd galaxy 
dSandage fc Tammannl I1981J 1, thus with essentially no 
bulge (see also Appendix [C]), but, from its broad nuclear 
emission li nes, has been observed to contain a Seyfert 1 
nucleus (jFilippeuko fc Sargen t UTosoh . A measurement of 
the mass of the central black hole in NGC 4395 helps our 
understanding of the co-evolution of galaxies and black 
holes in multiple ways - by adding to the statistics of 
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MBHs in low-mass galaxies and thus understanding the 
M-er relation and its scatter at low masses, by showing 
that MBHs can be present even in the complete absence 
of a bulge, and by exploring the demographics of MBHs 
in spiral galaxies. 

The mass of the black hole in NGC 4395 has been 
estimate d by a variety o f methods. Using reverberation 
mapping. iPeterson et al.l (1200511 estimated the mass to be 
M,= 3.6 x 10 5 Mg. However, more recently. lEdri et all 
(|2012f ) estimated the black hole mass using reverberation 
mapping to find M.= (4.9 ± 2.6) x 10 4 M 0 ; the differ¬ 
ent mass estimates are almost entirely due to modeling 
assumptions. The conversion of a time delay to a black 
hole mass is based on a virial factor dependent on the 
broad-line region (BLR) geometry. This virial factor is 
often calibrated by assuming that reverberation mapping 
black holes masses follow the same M.-er relation as those 
in inactive galaxies for which B H masses can be measured 
from stellar dynamics ('e.g. lOnken et al.ll2004f ). The esti¬ 
mates for NGC 4395’s black hole mass therefore rely on 
two assumptions: (1) that the M.-cr relation for quies¬ 
cent, mostly early-type galaxies is still valid for the spiral 
galaxies typically probed by re verberat ion m appin g, de¬ 
spite evidence of the contrary (|Greene et al.l 12010 ). and 
(2) that the BLR geometry remains similar at the lower 
mass inferred in NGC 4395. Evidence for a different 
mode of accretion has been sug gested by NGC 4395’s X- 


ray luminosity and variability 

Moran et al. 19991 200! 

but see Nardini & Risaliti 2011 



Obtaining a mass measurement of the MBH in NGC 
4395 in an independent and more direct way, such as dy¬ 
namical modeling, is therefore important. However, the 
masses of BHs in low-mass galaxies are much harder to 
measure than in high-mass galaxies, partly because of the 
lower surface brightness of the host galaxy and because 
the sphere of influence, in which the gravitional potential 
is dominated by the mass of the black hole, is smaller, 
but often also because of radially varying stellar mass- 
to-light ratios (M/L) due to recent star form ation or the 


presence of a nuclear star cluster (NSC. Walcher et al. 

2006 

: Rossa et al. 20061; Seth et al. 2006; Carson et al. 

2015 

Se 

mass 

ics 

)• 

veral authors have inferred black hole 
es from modeling the nuclear gas kinemat- 
in nearby massive galaxies (Macchetto et al. 

1997); van der Marel & van den Bosch 

1998c 

Verdoes Kleiin et al. 200(1 Neumaver et al. 

2007'; 

Barth et al.l 200lie iSarzi et al.l 2001; Cappellari et al. 

2002; iShaniro et al. 2006h. Here we nresent a dynamical 


measurement of the mass of the MBH and nuclear 
star cluster in NGC 4395, using high-spatial resolution 
adaptive optics assisted observations of the molecular 
gas kinematics from Gemini/NIFS. This is the l owes t 
mass MBH and, besides M60-UCD1 dSetlr et al.l 12014 b 
the lowest mass galaxy for which the black hole mass 
has been measured by dynamical modeling. 

The structure of this paper is as follows. In Sec¬ 
tion [21 we describe the photometric data, which we use 
to determine the stellar mass density (Sec. [4}, and 
the spectroscopic data, from which we determine the 
combined dark and luminous mass density by model¬ 
ing the gas kinematics (Sec. [5]). Throughout this pa¬ 
per, we assume that the distance to NGC 4395 is 4.4 
Mpc, which is in between the Cepheid distance (depend¬ 


ing on the choice of the per iod-luminosity relation ei¬ 
ther 4.0 or 4.3 Mpc, iThim et al.l I2004T ) and tip of the 
red giant branch (TRGB) distances (4.6 and 4.76 Mpc, 
IKarachentsev et all 120031 IJacobs et al.l 120091) . 

2. DATA 

2.1. Imaging data 

The center of NGC 4395 was imaged with HST/WFC3 
in 7 filter bands, ranging from the ultraviolet (UV) to 
near-infrared (NIR). Instead of a broad Vband, we used 
the medium-wide F547M filter, to exclude most bright 
emission lines. The observations are summarized in Ta¬ 
ble ffl Each filter band consisted of at least 4 exposures 
with different dither positions. These data are part of a 
larger survey to infer the formation mechanisms of the 
nearest nuc lear sta r clus ters a nd have already been pre¬ 
sented in iCarson et al.l (I2015I L however, to obtain opti¬ 
mal spatial resolution and better noise characteristics, 
we use a slightly different procedure to reduce the data. 

We drizzled the flat-fielded optical/UV images using 
the drizzlepac package in PYRAF. We found good re¬ 
sults in terms of resolution and noise by drizzling the 
optical images with a Lanczos3 kernel (with pixfrac set 
to 1 pixel). The World Coordinate System of the im¬ 
ages was aligned to that of the F438W band using the 
tweakreg package, after which the aligned flat-fielded im¬ 
ages were re-drizzled. The spatial alignment between the 
optical/UV bands is better than 0.1 pixel. 

Although we also imaged NGC 4395 in the NIR with 
HST/WFC3, we did not use these data for this work, 
since they have slightly lower resolution, whereas the sep¬ 
aration of the AGN contribution to the light of the NSC 
requires superior resolution. 

In addition to the WFC3 data, we also use archival 
HST/ACS F606W and F814W wide field camera (WFC) 
data. The F606W filter is important for identifying line 
emission, since it contains Ha. The F814W ACS data 
is necessary to make a F606W-F814W colormap. We 
drizzled the F606W and F814W data on the same grid 
as the WFC3 data, so that we were able to use these data 
to mask emission line regions in the WFC3 images. 

2.1.1. Point spread function 

Point sp read functions (PSFs) were generated using 
TinyTim (IKrist et al.l 1201 111 . This program generates 
raw PSFs, with the same geometric distortion as the 
flathelded HST/WFC3 images. Since we fit our Gal- 
fit models on the distortion-corrected drizzled images, 
we have to drizzle the PSF in the same way as the flat¬ 
helded data. For each band, we created therefore four 
intermediate PSFs with TinyTim , with the same pixel 
offsets as the dither pattern of the actual observations. 
These intermediate PSFs were injected in the hat-fielded 
science frames, which we filled with zeroes hrst. We driz¬ 
zled these frames with the same parameters as the other 
imaging data, making sure that no sky subtraction was 
performed during the drizzling process. We found that, 
after drizzling, the final PSFs were not centered on a 
pixel center, and we therefore shifted the four offsets that 
were used as input for TinyTim until the drizzled PSFs 
were centered on a pixel center. 

2.2. Spectroscopic data 
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Table 1 

Summary of HST/WFC3 observations 


Band 

Expose time 
[seconds] 

Zeropoint 

[mag] 

A\ 

[mag] 

F275W 

8 x 150 

22.6322 

0.094 

F336W 

4 x 254 

23.4836 

0.077 

F438W 

4 x 148 

24.9738 

0.062 

F547M 

4 x 120 

24.7477 

0.050 

F814W 

4 x 100 

24.6803 

0.026 


References. — (1) The extinc- 

tion values A\ were obtained from 
ISchlaflv fc Finkbeineil Q201U ) 


NGC 4395 was observed in the K band with Gem- 
ini/NIFS, a near-infrared imag e-slicing integral field 
spectrograph on Gemini-North ([McGregor et al.1 120031 ) 
with ALTAIR laser guide star adaptive optics. A nearby 
star was used to tune the slow focus for these obser¬ 
vations. Data was taken on three separate nights in 
2010: March 28th, May 5th, and May 22nd. Exposures 
were taken in a object-sky-object sequence and a total 
of 12 x 600s on source exposures were taken. A telluric 
calibrator with spectral type of A0V was taken on each 
n ight and used to remove telluric absorption. 

iFilippenko fe Hoi (|2003l ) tried to measure the velocity 
dispersion of the NSC to constrain the black hole mass, 
but could only obtain an upper limit on the dispersion of 
30 kms -1 . To get a better measurement of the disper¬ 
sion, we observed the NGC 4395 nucleus with NIRSPEC 
at the Keck II telescope. Unfortunately, the CO band 
heads, which are typical of evolved stellar populations, 
were absent, probably because of the continuum emission 
of the AGN, so that we were unable to obtain a disper¬ 
sion measurement from them. Since we do not use the 
NIRSPEC data for our analysis, we present the details 
of these data and the reduction procedures in Appendix 
E together with the details of our attempt to measure 
the dispersion from the CO bandheads in the NIFS data. 

2.2.1. NIFS data reduction and PSF determination 

The data was reduced as described by ISeth et abl 
(|2010D ; initial reduction was done using tools from the 
IRAF Gemini package that were modified to enable pro- 
pogation of the variance spectrum. The cubes were re¬ 
binned to a pixel size of 0'.'05 x 0'/05 using our own IDL 
version of NIFCUBE to make final data cubes for each 
night. The 12 data cubes were then combined after 
shifting the May data to match the March 28th data’s 
barycentric velocity. The line spread function (LSF) of 
the data was determined by reducing the sky cubes and 
combining them in the same fashion as the science data; 
the sky lines in each spatial pixel (spaxel) were then fit to 
Gaussians to obtain the variation in LSF across the chip; 
the full width half maximum (FWHM) of the LSF ranged 
from 4.0 A to 5.0 A across the chip with a median of 4.5 
A, a typical spectral resolving power of A/AA ~ 4800 
at the wavelength of the H 2 1-0 S(l) line used in our 
analysis. 

The intensity, velocity and velocity dispersion of the 
H 2 1-0 S(l), H 2 1-0 Q(l) and Bry lines were determined 
by fitting each line with a Gaussian. We found that close 
to the center of the cluster, the Bry line did not have a 
Gaussian shape, but could be fitted by two Gaussians: 


one with a narrow component (cr ss 40 kms -1 ) and one 
with a broad component (cr ~ 300 kms -1 ). The width 
of the broad component suggests the emission originates 
from the BLR. 

We used a Bayesian criterion to tell us where the data 
warranted a two-component fit to the Bry line. At these 
positions, the Bry line was fitted by a double Gaussian. 
To constrain the central potential of NGC 4395, we use 
the H 2 1-0 S(l), which has the highest S/N of the two 
H 2 lines. The H 2 1-0 S(l) kinematic maps are shown in 
Fig. [l] Although gas velocity dispersions can only be 
accurately determined for dispersions higher than ~15 
kms -1 , we note that our results are not sensitive to the 
dispersion of the disk (see §5). 

To model the PSF for the kinematics, we use the con¬ 
tinuum emission from the AGN. This is justified, since 
we know from the absence of the CO band heads that 
the continuum emission is made up almost completely of 
non-stellar emission (see Appendix lAl for more details). 
We have checked the width of the continuum against that 
of the flux from the broad line Bry emission and a star in 
the field. The FWHMs of these other components agree 
with the value of the cont inuum emis si on. 

We note that, as in ISeth et all (|2010h . the PSF 
can be modeled as a combination of a Gaussian (with 
FWHM=0"16), describing the instrumental PSF, and a 
Moffat profile, describing the seeing: £(r) = E 0 /[(l + 
(R/Rd) 2 ) 4 ' 765 )- We use Rd = 0"95. The inner core o f the 
P SF i s slightly broader than the one used by ISeth et abl 
(12010D . which is probably due to the nucleus of NGC 
4395, which was used for tip/tilt correction during the 
observations, being less bright. For the dynamical mod¬ 
eling, we decomp ose the profile into a multi-Gaussian 
expansion (MGE. lEmsellem et al.lfl994D . since it allows 
for an easy deprojection. 


3. STRUCTURAL PARAMETERS OF THE NSC 


The motions of the molecular gas in NGC 4395’s 
nucleus are determined by the gravitational potential, 
which depends on the mass of the MBH and the mass 
distribution of the stars around the MBH. There is no 
evidence that dense molecular gas in the center o f NGC 
4395 contributes much to the mass dBoker et all 120 111 ) . 
In order to derive the stellar mass distribution around 
the central black hole, we model the surface brightness of 
the NSC in the WFC3 images with a parametric mode l. 
For t his, we use the commonly used Sersic model (ISersid 
119681) . which is given by 


I(R) = I e exp {-&„ [(R/R e y - l] } 
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Figure 1. The H 2 1-0 S(l) velocity and velocity dispersion fields. North is up and East is left. White contours represent the intensity of 
the H 2 emission and black contours represent the total intensity in the Kg band. 


R e is the effective radius and I e is the surface brightness 
at that radius. The cuspiness of the profile is described 
by the Sersic index n. Wc will use the Sersic models of 
the NSC in the different filter bands to infer the stellar 
populations of the NSC (Sec. 0 ]), and to generate the 
stellar potential (Sec. 0 ]). To fit PSF convol ved Sersic 
mode ls to our data, we make use of Galfit ( iPeng et all 
l20f)l . 

The modeling of the NSC is complicated by several 
factors. First, part of the observed emission consists of 
continuum emission from the central accreting black hole. 
Because of the small spatial scale from which this emis¬ 
sion originates, we model this as a point source; the rela¬ 
tive flux in each band is however not well known. Second, 
a blue heart-shaped feature is clearly seen extending hori¬ 
zontally across the NSC in the F606W-F814W color map 
made from the ACS data (Fig. [2j. The similar F547M- 
F814W WFC3 data shows a much less dramatic color 
difference and, since the emission of this double-lobed 
structure is similar in morphology to the line emission in 
the NIFS data, we thus assume that this region is con¬ 
taminated by narrow line region emission from the AGN. 
Since the F606W filter contains [O III], Ha and H/3, this 
emission is more prominent in the F606W-F814W color 
map than in the F547M-F814W colormap. Although the 
emission line regions are brightest in our UV filters, they 
are clearly seen in the residual images of the other filters. 
There is thus no filter free of this emission. Third, th e 
NSC is quite compact (R c g ~ 4 pc; iCarson et al.1 2015), 
making distinguishing the AGN emission from that of 
the NSC somewhat degenerate. Of these complicating 
factors, the emission line region is the most difficult to 
deal with; in the following subsections we use three sepa¬ 
rate methods to try eliminating this component from our 
fits; comparison of these methods gives us some sense of 
our systematic uncertainties in the fits. We will then use 
these models in our black hole mass determinations in 
Section 0 ] 

For comparison, we also pro vide t he st ructural param¬ 
eters of the NSC obtained by ICarson et al.i ( 201 5) from 


the same WFC3 data, but without dealing with the emis¬ 
sion line region in a special way, in Table [2] The param- 
eters of the NSC as determined by ICarson et alJ ( 2015 ) 
are sufficient for analyzing the global shape of the NSC. 
However, since the black hole mass is only a fraction of 
the NSC mass, a slight alteration of the density profile 
of the NSC by the emission line region could change the 
amount of mass at the center of the cluster, and thus 
change our black hole mass determination. A very accu¬ 
rate treatment of the emission line region is thus crucial 
for our dynamical modeling. 

3.1. Masking emission-line regions 

The most obvious way to get rid of the emission line 
regions is by masking them out. We identify all central 
pixels in the F606W-F814W color map that are bluer 
than 0.38 mag, and create a bad pixel map for Galfit 
based on these pixels. 

We fit the light distribution of the NSC with a com¬ 
bination of a Sersic profile and a point source. We first 
determine the Sersic index and effective radius in the 
F547M band, and then keep those fixed between differ¬ 
ent bands. We find a Sersic index n = 2.25 and R e g ~ 3.6 
pc. We note that these values for the Sersic are not atyp¬ 
ical for NSC s: the Milky Way NSC h as n ft 3 and M32 
has n ft 2.3 ((Graham & Suitie'ill2009( ). 

3.2. Iterative residual subtraction 

Instead of masking the regions, we can also try to 
model and subtract the emission line contribution in each 
filter. We tried two ways to do this. The first method, 
which we dub “iterative residual subtraction”, makes use 
of the fact that the emission line regions show up as 
strongly positive residuals after fitting. For each band, 
we identified strongly positive residuals from our Gal¬ 
fit fit, with the bluest pixels still masked, as emission 
line regions. The thus created “emission line map” was 
then subtracted from the original image, which was then 
fit again with Galfit. This procedure was iterated a 
few times. It is a priori not clear that this procedure 
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only identifies the non-axisymmetric emission, since fit¬ 
ting with the wrong NSC model can lead to the removal 
of stellar emission too. We were therefore very careful 
in identifying at which intensity to clip, which unfortu¬ 
nately leads to some arbitrariness. 

The final structural parameter data show a somewhat 
more compact NSC, with Sersic index n = 4.4 and effec¬ 
tive radius R e s = 3 pc. The NSC is 0.07 mag brighter 
in the F814W band than with the masking method. Al¬ 
though it may seem counter-intuitive that subtraction of 
the residuals would lead to a brighter NSC, the overall 
curvature and compactness of the surface brightness pro¬ 
file are significantly changed, resulting in this change in 
brightness. 

3.3. Subtraction of the emission-line regions 

Inspection of the F275W image suggests that the nar¬ 
row line emission contamination is most prominent in the 
F275W band, possibly due to the lack of stellar emission 
in these bands. Its isophotes are the most non-circular, 
and resemble those of the Bry emission in the NIR. In¬ 
stead of masking the emission line regions, we therefore 
attempt to remove the emission line regions by subtract¬ 
ing the F275W image from the redder bands after ap¬ 
propriate scaling. Since the diffraction pattern in the 
F275W band is different from the other bands, we first 
subtract the central point source from the F275W band. 

To scale the emission of the F275W band to the other 
bands, we assume that the emission line region is the 
narrow line emission of the AGN. We therefore use a typ¬ 
ical narrow line emission spectrum, from t he we ll s tudie d 
nearby Seyfert 2 galaxy NGC 1068 (jSpinelli et al.ll2006h . 
to determine the flux ratio between the F275W band and 
the redder bands. Subtracting the scaled point source 
subtracted F275W image from the other bands gave sur¬ 
prisingly good results in the F547M band whereas the 
F814W band required another 5% scaling in the flux of 
the F275W band. 

The best fitting Sersic profile, determined in the 
F547M band, has Sersic index n = 1.4 and i? e ff = 4.4 
pc. The luminosity of the NSC is lower than for the 
other two methods: compared to the masking method 
the NSC is 30% fainter in the reddest band and in the 
bluer bands even more. We note that it is possible that 
some of the subtracted F275W emission is in fact stellar. 

3.4. Preferred structural model 

We summarize the structural parameters found for the 
four different colors using the t hree differ ent methods, to¬ 
gether with the results of [Carson et al.l 1 2015 1. in Table 
[2 To assess which method describes the cluster profile 
the best, we compare how the models fit the outer, un¬ 
masked portion of the NSC. We avoid consideration of 
the inner parts of the fit, because the iterative residual 
subtraction gives a low y 2 by construction. We therefore 
decided to calculate a y 2 , with 

2 _ (I» — Models) 2 

X ~ 2^ 0-2 Cl 

i£Good pixels 1 

where the good pixels are all unmasked pixels (deter¬ 
mined from the bad pixel map of the masking method) 
in an annulus of 4-20 pixels from the cluster center, cor¬ 
responding to roughly 1 5 Here I,;, Models and cy 


denote the observed intensity, model intensity and obser¬ 
vational uncertainty of each pixel value. 

The masking and iterative residual subtraction method 
have the lowest y 2 per pixel values (4.3 and 4.1), the y 2 of 
the F275W subtraction method is slightly higher (4.7). 
Given the small difference in performance between the 
masking method and the iterative subtraction method, 
we adopt the simpler approach and use the results from 
the masking method as our preferred NSC parameters in 
our dynamical modeling. To get a sense of the system¬ 
atic uncertainty in our BH mass measurements due to the 
poorly constrained NSC profile, we also perform dynam¬ 
ical modeling using NSC parameters from the other two 
methods. The profiles of all uncolvolved Sersic models of 
the NSC are shown in Fig. (3[ as well as a PSF-convolved 
model including the central point source emission. 

4. STELLAR POPULATION SYNTHESIS MODELING 

An important constraint on the mass of the nuclear 
star cluster comes from the analysis of its stellar pop¬ 
ulations. Although the magnitudes of individual bands 
vary depending on the fitting method and masking, some 
of the colors are surprisingly constant. In particular, 
the F547M-F814W color is 1.0, 0.87, and 0.96 for the 
three methods that do deal with the emission line re¬ 
gions. The F336W-F438W colors show however a bigger 
spread with F336W-F438W = -0.21, -0.19, and 0.33. 

To constrain the mass of the NSC, we need to infer its 
stellar mass-to-light ratio (M/L) in a certain band. A 
common choice for this is the reddest band, since it is 
least sensitive to the presence of young stellar popula¬ 
tions and extinction. We therefore infer the M/L ratio 
of the cluster in the F814W band. 

We make use o f the Flexi ble Stellar Population Synthe¬ 
sis code CFSPS. [Conroy et al.l l2009h to calculate colors 
and masses of stellar populations with a range of metal- 
licities and ages. We select four metallicities, which cover 
a very metal poor scenario (log Z/Zq = —1.50), a metal 
poor scenario (log ZfZ & = —0.69), solar metallicity 
(log Z/Zq = 0.0) and super solar metallicity (log Z/Zq = 
+0.20), and 180 different ages, logarithmically spaced be¬ 
tween 100 Myr and 15 Gyr. 

We assume the initial mass function (IMF) of [Kroupai 
(|2001f ). FSPS takes into account that the mass of the 
cluster evolves over time; we assume that all stellar rem¬ 
nants are contained within the cluster, but do not calcu¬ 
late the mass in gas, since it is not clear that this gas is 
still present. We then let FSPS make photometric pre¬ 
dictions for four of our five optical bands. We exclude the 
F275W because of the large fraction of likely non-stellar 
emission. 

4.1. Single stellar population (SSP) fits 

As a first attempt, we fit an SSP to the magnitudes 
derived with the Galfit fits. The photometric uncer¬ 
tainties on the data are small, typically 1 per cent. This 
uncertainty does not reflect the systematic errors such 
as model mismatch, or remaining line emission. The un¬ 
certainty represent therefore a lower limit on the real 
uncertainty. 

We use the y 2 as our figure of merit for the fits. For 
each age and metallicity in the library, we calculate the 
predicted magnitudes, given a stellar mass. A mass-to- 
light ratio is then calculated by dividing the F814W lu- 
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Figure 2. Color maps of the nucleus of NGC 4395. On the left, we show the F547M-F814W color map based on HST/WFC3 data. On 
the right, we show the F606W-F814W color map based on HST/ACS data. The asymmetric blue emission is probably due to line emission. 


Table 2 

Best fit parameters 


Band 

Magnitude 1 

[mag] 

Effective radius 
["] 

Sersic index 

Point source magnitude 
[mag] 

Carson et al. 2015 

F336W 

17.61 

0.130 

1.86 

17.25 

F438W 

18.44 

0.168 

1.83 

18.22 

F547M 

17.72 

0.176 

2.80 

18.15 

F814W 

16.98 

0.221 

1.41 

17.03 

Emission line masking 

F336W 

18.42 

0.150 

2.25 

16.80 

F438W 

18.63 

0.150 

2.25 

17.93 

F547M 

17.90 

0.150 

2.25 

17.99 

F814W 

16.89 

0.150 

2.25 

17.10 

Iterative residual subtraction 

F336W 

18.17 

0.123 

4.43 

17.01 

F438W 

18.36 

0.123 

4.43 

18.23 

F547M 

17.69 

0.123 

4.43 

18.36 

F814W 

16.82 

0.123 

4.43 

17.32 

F275W subtraction 

F336W 

19.75 

0.183 

1.41 

17.55 

F438W 

19.32 

0.183 

1.41 

18.78 

F547M 

18.19 

0.183 

1.41 

18.66 

F814W 

17.23 

0.183 

1.41 

18.38 


Note. — (1) Magnitude of the Sersic component in the Vega system 


minosity by the mass, which can be different by a few 
per cent from the theoretical M/L ratio from the stellar 
library, since it allows for a statistical error in the F814W 
band luminosity. 

Fig. |7] shows a probability distribution function of the 
fitted M/Lp 8 i 4 w- The best fitting single stellar popula¬ 
tion is old and metal poor. The x' 2 of individual fits is 
not perfect, in the sense that the reduced xt > 1- An 
old metal poor population provides a better fit than a 
younger metal rich population. The age and metallic- 
ity of the nuclear star cluster in this scenario would be 
consistent with t he NSC being built up from ac creted 
globular clusters llLotz et al.ll2001fc lAntoninill2013lh 

4.2. Composite stellar population (CSP) fits 

Although it is possible that the NSC formed from old 
metal-poor globular clusters, sitting at the galaxy cen¬ 
ter it is very likely to have been e xposed to gas accre¬ 
tion and subsequ ent star formation dWalcher et al.ll2006b 
iSeth et all I2006T ). A more physical mass model of the 


NSC should therefore include the possibility of an addi¬ 
tional younger population. We therefore consider CSP 
fits where we assume that the NSC is built up from two 
components: an old and metal poor component, which 
probably contains most of the mass, and a younger com¬ 
ponent with equal or higher metallicity. For each possible 
combination of metallicities and ages, we calculate the 
best fitting weights of the two stellar populations from 
our four photometric bands and use this to obtain for 
each metallicity and age combination a likelihood and 
a M/L ratio in the F814W band. We show the nor¬ 
malized likelihoods of the fitted M/Lfsi 4 w ratios from 
composite stellar populations as a thick solid line in Fig. 
01 together with the 68% confidence interval. For SSP 
fits, M/Lf814W ratios up to 2 were allowed. For compos¬ 
ite populations, higher M/L ratios are also allowed. The 
reason for this is that it is now possible to fit the red part 
with an old metal-rich population, which has a slightly 
higher mass-to-light ratio than the old metal-poor pop- 
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Radius farcsecl 

Figure 3. Unconvolved Sersic models of the nuclear star cluster, 
derived from fitting the F547M band by masking the line emission, 
iteratively subtracting residuals and subtracting the scaled F275W 
image. We also show the data overplotted on the PSF-convolved 
Sersic + point source model that was derived by masking the line 
emission. The red vertical line shows the effective radius of this 
Sersic model of the NSC. 

ulation. We note that it is possible that the emission 
line contamination in the blue bands mimics a young 
stellar population. Whether this contamination is from 
young stars or emission lines is in reality not very rele¬ 
vant though, because neither would contain much mass. 

The best fit CSP models have a main component of 
solar metallicity and intermediate age (3 Gyr) combined 
with a young component. This would point at a scenario 
in which the nuclear stellar cluster was built up locally 
from enriched gas and continues forming stars, albeit at 
a low level. 

4.3. Empirical M/L ratios 

Observations of globular clusters suggest that SSP 
models with a standard IM F may not accurately pre - 
dict the M/L of old systems (IStrader et al.l 1200911201 111 . 
Although we do not know how the nucleus of NGC 4395 
formed, the colors of the cluster are consistent with those 
of globular clusters, thus motivating us to examine em¬ 
pirical constraints on the M/L from globul ar clusters. 

We use the sample of GCs in M31 of IStrader et al.1 
(|201 If ) for which dynamical masses have been calculated. 
We then look up the corresponding J-band ma g nitude s 
in the Revised Bologna Catalogue ([Galleti et al.l 120041 1. 
correct these for extinction, as in Strader et al., and con¬ 
vert the /-band magnitude to F814W magnitude. This 
yields 177 unique sources, for which we can determine 
the M/L in the F814W band. As in Strader et al., we 
use the virial mass as an indicator of the dynamical mass 
of the cluster. 

It is known that the mass-to-light ratio of globular clus¬ 
ters varies as a function of luminosity; low-mass globu- 
lars are sensitive to dynamical processes that can alter 
their M/L ratio dKruijssen fc Mieskel 12009 ). We there¬ 
fore chose to use only globular clusters heavier than 10 6 
Mq and with color 0.8 < V — I < 1.0, since our V-I 
color determination of NGC 4395’s cluster is the most 
robust and least affected by possible additional young 
stellar populations. We show the distribution of empiri¬ 
cal M/Lp 8 i 4 w ratios inferred for these 14 clusters in Fig. 



Figure 4. Normalized probabilities of the mass-to-light ratios 
of single and composite stellar populations. An SSP fit to the 
colors of our favorite structural model (thin black line) prefers a 
relatively metal poor and old population, because old metal-rich 
populations underestimate the emission in the bluest bands. The 
range of M/Lp 8 i 4 W ratios for composite populations (thick line) is 
broader, since metal-rich populations in combination with young 
populations can now also fit the data. Although scaled to a com¬ 
mon value for this plot, the likelihood of the CSP fit is higher than 
that of the SSP fit (Ax 2 ~ 7). The CSP fits to the the colors 
of the other two structural models are shown as a green dotted 
line (for the iterative residual subtraction method) and red dash- 
dotted line (for the F275W subtraction method). We also show as 
a histogram the empirical mass-to-light ratios of M31 GCs in the 
I band, based on data from ^trader et al.l H2011J) and the Revised 
Bologna Catalogue 

HI 

The empirical M/Lfsi 4 w ratios of the GCs vary from 
0.5-3.0 Mq/Lq and peak around between 1-1.5 M 0 /L 0 , 
consistent with the M/Ls from the theoretical stellar 
populations analysis. We note that including GCs with 
colors outside our chosen color range does not change 
the peak value or width of the histogram. As the sources 
shown here are the most massive GCs, it is possible that 
some of thes e contain MBHs and therefore have e levated 
M/L ratios dMieske et al.1 120131 : iSeth et al.l 120141) . Nev¬ 
ertheless, we can conclude that the empirically derived / 
band M/L ratios of massive GCs and the M/L ratio of 
NGC 4395’s NSC derived from stellar population synthe¬ 
sis modeling are in excellent agreement with each. There¬ 
fore, for our dynamical modeling we use the constraint 
from the CSP modeling that the M/Lfsi 4 w is 1-3 ± 0.6 
(lcr uncertainty). 

5. GAS DYNAMICAL MODELING 

The left panel of Fig. Q] shows the line-of-sight velocity 
of the molecular hydrogen gas. The kinematics show a 
clear rotational signal, with a position angle (P-A.) of 
~ —75° that is almost orthogonal to the P.A. of ~ 28° 
found for the radio jet dWrobel fe Holl2006lf . 

To find the best fitting potential, we adopt the com¬ 
monly used assumption that the gas is on c i rcular orbits 
around the black hole dBertola et al.l 09981: iBarth et alJ 
120011 : iNeumaver et al.1 l2007t ISeth et al.l 120101) . In this 
case, the circular velocity of the gas is only proportional 
to the enclosed mass. We assume a spherically symmet¬ 
ric density profile, but note that the NSC appears to be 
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Table 3 

Multi-Gaussian Expansion of the light profile of the 
NSC in the F814W band. 


Component number 

Surface luminosity 
[L 0 /p C 2 ] 

Dispersion 

["] 

1 

92502 

0.00082 

2 

102553 

0.0027 

3 

94608 

0.0076 

4 

70342 

0.019 

5 

40710 

0.042 

6 

17737 

0.086 

7 

5653 

0.16 

8 

1288 

0.28 

9 

206 

0.47 

10 

23 

0.75 

11 

1.7 

1.14 

12 

0.08 

1.71 

13 

0.001 

2.65 


slightly flattened (q « 0.95). We characterize the light 
profile of the cluster usin g a s pherical Multi-Gaussian 
Expanson (MGE) model (iBendinellil 11991 ). We create 
the MGE using the code from iCappellari (1200211 to fit 
the 1-D surface brightness profile of our Galfit models 
in the F814W band. The MGE of our prefered structural 
model is tabulated in Table [3] 

We then derive the stellar density profile, p(r ), by de- 
projecting the MGE of the stellar surface brightness pro¬ 
file in the F814W band, E(i?), and multiplying by a dy¬ 
namical M/L ratio: 


N 


w) = E 


L, 


i= 1 


27rcr- 


■ exp 


-R 2 


2af 


N 


»M=£f^exp 
i= 1 ( 27r ) 2 W 


„ 21 


2 ai 


( 2 ) 

( 3 ) 


where R and r are the projected and intrinsic radius, 
Li the luminosity of each Gaussian, (M/L) the dynam¬ 
ical mass-to-light ratio of each Gaussian (which varies 
between dynamical models but is the same for all Gaus- 
sians of the MGE) and cq the dispersion of the Gaussian 
component. 

The gas disk orbiting in this potential is built up from 
rings, each with its own position angle and major axis, 
but all concentric. We assume a single inclination for 
all rings (see below). Neither the kinematic observations 
nor the H 2 gas morphology provide any indication that 
the disk is eccentric. __ 

We use Kinemetry (|Krainovic et al.ll2006f) to analyse 
our two-dimensional molecular gas velocity field. The al¬ 
gorithm assumes the velocity field can be described by a 
series of ellipses with varying maximum line-of-sight ve¬ 
locity, position angle (P.A.), cllipticity, and semi-major 
axis. Along each ellipse, the velocity is modulated by 
a cosine term. Kinemetry determines the best fitting 
set of concentric ellipses and provides us thus with a line 
of nodes (semi-major axis, PA, q) and the maximum ve¬ 
locity of each ellipse. The line of nodes is used as an 
ansatz in our modeling code to generate a model veloc¬ 
ity field. During the course of this work, we found that 
our models performed better with a fixed axis ratio for all 
ellipses. We therefore use a single inclination for all rings. 
For each radius, we calculate the circular velocity using 


v\ = rd<&/dr. For the MGE, this equation becomes: 

v 2 _ ^, -Mbh + ATmge (r) 
r 

with Mmge(i') given by Eq. 49 in ICappellaril (120081 !. 
Given a ring’s inclination and position angle, the line- 
of-sight velocity at each position along the ring is then 
calculated by multiplying the velocity with the cosine of 
the angle along the ellipse (the true anomaly) and the 
sine of the inclination. We determine the values of the 
line-of-sight velocity for a 10 x oversampled grid in (x, y) 
coordinates by interpolating between the rings. 

Since the H 2 kinematics are well described by single 
Gaussians, with dispersion values below the resolution 
of the spectrograph, we decided to predict the intensity 
weigted mean velocity ( Iv) instead of the peak velocity of 
a line with a possibly complex structure. Besides a model 
velocity field v m (x, y), each model also has an associated 
broadening field with it, cr m (x,y), which describes the 
intrinsic dispersion of the H 2 gas. Inside the telescope, 
the light distribution of each frequency is convolved with 
the PSF. The observed data cube Z can thus be described 
by 

Z(x,y,v) = JJ I(x,y)(j>(x,y,v) PSF(x' - x,y' - y) dx'dy 

( 5 ) 

with I(x, y) the (deconvolved) intensity field of the source 
and <j>(x, y, v) a position dependent function, 


<t>(x,y,v) = 


\j2iT<r m (x,y) 


exp 


(v - Vm(x,y )) 2 


2 <Tm( X >y) 


( 6 ) 


describing the velocity v of the disk at location (x. y) and 
the thermal and turbulent broadening of the velocity field 
a m . The mean velocity can thus be found by integrating 
Z(x,y,v) ■ v. Since v is not dependent on the position, 
we can change the order of integration: 


(Iv) = J Z(x,y,v)vdv 

= JJ I(x,y)Vm(x,y)PSF(x' - x,y' - y) dx'dy'. 

( 7 ) 

To simulate the NIFS observations, we can thus mul¬ 
tiply the velocity field by the oversampled deconvolved 
intensity field of the H 2 gas, and convolve this with the 
oversampled NIFS PSF. We rebin the convolved field to 
the size of NIFS pixels and divide it by the convolved in¬ 
tensity field to find the predicted velocity for each NIFS 
pixel. 

We note that this is computationally a significantly 
faster way to calculate the first moment than the full line 
modeling, since this involves doing several PSF convolu¬ 
tions for all velocities. However, this faster method will 
break down when the velocity at the peak of the emission 
line (which is approximately what we measure by fitting 
a Gaussian during the NIFS data analysis) is no longer 
close to the average velocity of the line (as defined by 
Eq. [71). 

If the gas disk is partly supported by internal pressure, 
the black hole mass in the cold disk scenario is underes- 
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Figure 5. Maximum line-of-sight velocity (u ro ^ sin(z)) and ob¬ 
served velocity dispersion (a) as a function of radius, determined 
with Kinemetry. Values of a below 15 are generally not reliable, 
as these values are well below the instrumental resolution of NIFS 
(FWHM = 60 kms -1 ). Also shown is the velocity corrected for 
the best fit inclination (^ ro t)- 

timated. From Fig. [5j it appears that pressure could 
be contributing significantly. However, the instrumental 
dispersion of NIFS is cri nst ~ 25 kms -1 ; from previous 
measurements we have found that dispersion values be¬ 
low 15 kms -1 are not reliable. The non-central areas 
with the highest S/N show also the lowest dispersion val¬ 
ues, close to 5 kms -1 , indicating that the dispersion in 
the other parts may also be lower. Full line modeling 
of our best fit model shows that there is no need for an 
additional pressure term. 

5.1. Fitting 

Our dynamical models have the following free param¬ 
eters: 

1. The dynamical M/L of the nuclear star cluster in 
the F814W band. We use a logarithmic scale be¬ 
tween 0.5 and 5. 

2. The mass of the IMBH. We assume that this mass 
is between 10 3 and 10 7 M 0 logarithmically scaled. 

3. The inclination of the system. Although Kineme¬ 
try provides us with an ellipticity of each ring, 
which we convert to an inclination, we also try out 
configurations with a constant (possibly negative) 
angle added to all rings, and with all rings fixed to 
a single inclination value. 

Each point in this parameter space corresponds to a dy¬ 
namical model. The likelihood of each data pixel is cal¬ 
culated from the squared difference of the observed and 
predicted velocity divided by the uncertainty of the ob¬ 
served velocity. In Fig. [I] a bright region is visible west of 
the nucleus, at offset 0.5-1.0 arcsec. We exlude therefore 
all pixels with X offset higher than 0.5 arcsec from the 
fit. We also exclude velocities with uncertainties higher 
than 5 kms -1 . For these uncertainties, fluctuations in 
the sky can cause systematic errors in our velocity de¬ 
terminations. This mainly affects isolated pixels of low 
S/N in Fig. m 


5.2. Results 

We created dynamical models for 25 different inclina¬ 
tions, each sampling 25x25 models in black hole mass 
and M/L. The M/L and black hole mass were sampled 
logarithmically. Fig. |T]shows the likelihood of these mod¬ 
els, projected by taking the maximum likelihood along 
the inclination axis or the M/L axis. Our dynamical 
models seem to prefer a configuration in which the gas 
disk has M/Lpsuw = 1.3 and inclination of 37°. The 
M/Lp 8 i 4 w value is in good agreement with the stellar 
population analysis above. The reduced y 2 of our best 
fit model is 3.9. A y 2 value > 1 is not unexpected 
with our simplifying modeling assumptions. We there¬ 
fore scale the y 2 of our models so that the best fit model 
has y 2 = 1- For our 3 a confidence limits, we select 
models for which the rescaled y 2 is within Ay 2 = 14 (3 
degrees of freedom) of our best fit model’s y 2 . 

We have combined the results of the composite stellar 
population models with dynamical constraints. In Fig. 
[HI we show the results of this combination. The 1-cr con¬ 
tours do not change much compared to the fits based on 
dynamical models only, since the best fit dynamical mod¬ 
els agreed almost exactly with the best fit stellar popula¬ 
tion synthesis models, the joint analysis helps to exclude 
models with low inclination and low black hole mass. 
From the combined modeling of the stellar populations 
of the NSC and the dynamics of the molecular gas, we 
find a best fit black hole mass of M = 4jj® x 10 5 M 0 (3cr 
uncertainties). The velocity field of the best fit model is 
shown in Fig. [ 6 j For the parameters of this model, we 
also did the full line modeling, so that we could calculate 
the dispersion field. The central dispersion is accurately 
reproduced by this model, together with the depressions 
in dispersion along the H 2 disk, without adding any in¬ 
trinsic dispersion for the gas. This shows that our mod¬ 
eling of a cold disk is justified. 

The modeling results for different structural parame¬ 
ters of the NSC are shown in Figs. [Til and fl2l in the Ap¬ 
pendix. For the iterative residual subtraction method, we 
find a BH mass similar of M = 3+2 4 x 10 s M 0 , which is 
very close to the mass found using the masking method. 
The F275W subtraction method for determining the po¬ 
tential gives a higher black hole mass of M = 7+® x 10 s 
M 0 . 

6. DISCUSSION AND SUMMARY 

6 . 1 . Uncertainties and systematic errors in our mass 
determination. 

Several components of our modeling procedure are 
prone to systematic errors; the structural parameters of 
the cluster are uncertain, we are make simplifying as¬ 
sumptions about the shape of the potential, and it is 
possible that the observed gas does not be trace the grav¬ 
itational potential. 

The largest uncertainty comes from our poor knowl¬ 
edge of the structural parameters of the cluster. Despite 
this, dynamical modeling with the different structural 
models from Sec. [3] for the potential of the NSC yields 
best fit black hole masses that agree with our preferred 
model to within 1-cr. However, our analysis relies criti¬ 
cally on the assumption that we can separate the contin¬ 
uum emission of the AGN, which is Upsuw ~ 1-3 x 10 6 
L 0 from the stellar emission of the cluster Lp’ 8 i 4 iY ~ 
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Figure 6. The same as Fig. \T\ but now for the model velocity and dispersion field. In the bottom panels we show a horizontal cut 
through the velocity and dispersion fields. The data are shown with black error bars, and the model fields are shown with a red dashed 
line. Although the dispersion was not fit in our models, the peak is reproduced at the right magnitude. Note that we consider values for 
the velocity dispersion a < 15 kms _1 to be unreliable. 


1.6 x 10 6 L 0 . 

The kinematic signature observed in the H 2 gas does 
not necessarily have to be a rotating disk; it could also 
be tracing a weak outflow. However, there are a number 
of lines of evidence that suggest this is unlikely. 

1. The H 2 emission in other Seyfert gal axies is known 
to trace gas disk components ('e.g. lRiffel et al.ll2009i : 


IHicks et a!1l2?Ml20l3lt . 

2. The analysis of the velocity field with Kinemetry 
shows that the velocity along each ellipse follows 
almost a perfect cosine. An outflow would not nec¬ 
essarily have the same signal. 

3. IWrobel fe Hoi (|2006f present radio continuum ob- 
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Figure 7. Results of dynamical models while using photometric model B. The plots show the likelihood of the data for models with 
free inclination, M/L ratio and black hole mass plotted Left: as a function of M/L ratio and black hole mass, maximized for inclination. 
Right : as a function of inclination and black hole mass, maximized for M/L. Contours encompass models within 1 and 3cr distance from 
the best-fit model. 


Masking method 



L°g(M BH [M 0 ]) 

Figure 8. Results of dynamical models combined with constraints 
from stellar population synthesis models. 


servations of the central parsecs of NGC 4395 and 
find an elongated structure, possibly an outflow, 
along a PA of 28 degrees, more or less orthogonal 
to the rotation axis of the H 2 disk. The orientation 
of this jet structure suggests a disk interpretation 
for the H 2 gas instead of an outflow. 

In our models, we assume that the cluster is spherical 
and that the gas is not pressure supported. Loosening 
the first assumption could lower the black hole mass by 
a few percent. Introducing a pressure component will 
increase the black hole mass. 


6 .2. Comparison to previous measurements 

Wit h the exception of M3 2 (jVerolme et a,l .1120021) . NGC 
7457 (iGebhardt et al.ll2003H and M60-UCD1 ( Seth et all 
l20ll . very few dynamical mass measurements of massive 
bla ck holes in galaxies with low velocity dispersion s ex¬ 
ist (|Kormendv fa Holl2013t iMcConnell fa Mall2013f) . For 
most galaxies with dispersion below 100 kms -1 , masses 


have been determined using the reverberation mapping 
method. The formula to determine the mass is usually 
written as: 


Mbh = f 


(cAt)(Av) 2 

G 


( 8 ) 


where At is a delay time between variation of continuum 
and line. The first bracketed term is thus a measure of 
distance from the black hole. Av a proxy of the average 
rms velocity of the line-emitting gas. For non-Gaussian 
line shapes, it is not clear what quantity should be used 
here, FWHM or a. Virialization normalization factor / is 
usually called the f-factor, and is often calibra ted against 
t he M -cr relation for high-mass AGN (but see iHo fa Kiml 

l20ll . 

Three mass estimates for the black hole in NGC 4395 
were carried out using this technique, but with different 
data. iPeterson et al.: ( 2005 (1 used the C IV line in the UV 
to estimate the time lag, and converted this to a black 
hole mass using / = 5.5 (based on high-mass AGN) and 
using the velocity dispersion <7 of the line for Av (which 
in this case was very similar to the FWHM). Th i s mass 
measurement was confirmed by iDesroches et all (1200611 
using Ha RM, although their measurement was not sig¬ 
nificant. lEdri et gill ( 2012 ) used broad-band reverbera¬ 
tion mapping, which targets an individual line, in this 
case Ha, but uses a broad band filter to measure the 
variation in the line. These authors used / = 0.75, based 
on a simple virial estimate of a spherical BLR, and used 
instead the FWHM of the line. Our black hole mass es¬ 
timate is in agreement with the estimate by Peterson et 
al, suggesting that a high f-factor may be applicable in 
this particular case. 

In addition to the RM estimates, several additional 
BH mass estimates have been made for NGC 4395. Sev¬ 
eral of these studies use estimates of the BLR size from 
the AGN continuum, and then use a line width to de¬ 
rive a mass. iLira et akl ( 199 9) estimate the black hole 
mass to be 5 x 10 5 M 0 (assuming a distance to NGC 
4395 of 5 Mpc), based on the width of the H/3 line and 
their estim ate of the bolometric luminosity of the AGN; 
iKraemer et al.1 (1199911 estimate a mass of 1.5 x 10 5 M 0 (at 
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2.6 Mpc distance) based on H/3 and a size estimate of 
the BLR from photoionization models. Using H/3 an d 
the BLR size-luminosit y relation of lKaspi et al.l (2000), 
iFilinoenko fo Hoi (l2003l l estimate a mass 10 4 M@. 

A relation between black hole mass, bolometric lumi¬ 
nosity and break frequen cy of the X-ray power spectrum 
has also been developed dMcHardv et al.ll2006lf . In NGC 
439 5 the break frequency was measured to be ~ 2 x 10~ 3 
Hz (IVaughan et al.ll20051f. If we use the bolometric lumi¬ 
nosity from lLira et all (119991 (1.2 xlO 41 erg s' 1 ), we find 
a black hole with mass M = 5 x 10 4 M 0 , which is sig¬ 
nificantly lower than our estimate. However, there are 
some indications that the X-ray luminosity is unde res- 
timated due to extinction (iNardini fe Risalitl 1201 111 , so 
that it is possible that the bolometric luminosity is in 
fact higher, which would increase NGC 4395’s black hole 
mass estimate from the McHardy relation. 

6.3. A well-measured black hole mass in a bulgeless 

galaxy 

NGC 4395 is currently the only bulgeless galaxy with a 
dynamically measured black hole mass. Here we compare 
this galaxy to two similar bulgeless galaxies that appear 
to lack MBHs. 

At a distance of 4.4 Mpc, NGC 4395’s absolute mag¬ 
nitude is Mb = —17.65 mag and My = —18.10 mag, 
slightly less luminous than M33, and similar to NGC 
4244, an highly inclined Sc galaxy which lies at the 
same distance from Earth. NGC 4244 and NGC 4395 
have si milar maximum rotation velocitie s ( V max ~ 90 
km s' 1 , iSwaters et al.l fl999t lOllind 179961) which are 20 
km s' 1 lower than M33’s. Although all three galax¬ 
ies have similar morphological properties, NGC 4395 
is the only galaxy with an obvious MBH. The up¬ 
per limit on the mass of the MB H in NGC 4244 is 
5 x 10 5 M 0 (|De Lorenzi et al.1 120131 ) and for M33 is as 
low as a few thousand solar masses (Merritt et al.ll2001t 
iGebhardt et al.l 1200lH . The_ masses of t he NS Cs in 
M33 (1.4 — 2.1 x 10 6 M 0 , iKormendv fc McClurel 119931 : 
lHartmann et al.ll2011fl and NGC 4395 are essentially the 
same, wh ereas the NSC of NGC 4244 is 5 times heavie r 
(1.1 x 10 7 lHartmann et al JI201 1( : iDe Lorenzi et aDl2013ll . 
Despite having a black hole at its center, NGC 4395 has 
succeeded in producing an ordinary NSC with a_ mass 
typical of those in late type galaxies (iSeth et aJI 12008T) . 
NGC 4395’s /-band central surface brightness is, as ex¬ 
pected, slightly lower than M33’s. The only obvious dif¬ 
ference between NGC 4395 and these other two galaxies 
is the kinematic lopsidedness of NGC 4395. Lopsidedness 
is however not an uncommon featur e among late-type 
spiral galaxies dMatthews et al.lfl998T) . 

It is interesting to see how NGC 4395’s black hole 
mass fits in with scaling relations of MBHs in more mas¬ 
sive galaxies. NGC 4395 is not a dispersion dominated 
galaxy, and we do not know the value of its velocity dis¬ 
persion. Comparing it to the M-cr relation would t hus not 
be very useful. However, IKormendv fe Hoi (12013T) derive 
a relation between black hole mass and the A-band lu¬ 
minosity of the bulge (M-L), with which we can compare 
our measurement. In Appendix ICl we perform a photo¬ 
metric decomposition of NGC 4395. Although we believe 
that the central upturn in the light profile of NGC 4395 
is due to a bar, we can use the luminosity of the bar as an 
upper limit on the luminosity of the bulge, and compare 


the black hole mass with this upper limit. We do not 
know the A-band luminosity of the bar. To convert our 
/-band measur ement of L = 4.1 x 10' Lq, to A band, we 
use i — I = 0.4 (iJordi et al.ll2006lf and I—K = 2, which is 
common for old stellar populations, and is p r obab ly red¬ 
der t han the actual c olor of the bar (IPeletier fe de Griia 
I1998T) . Formula 2 in IKormendv fe HoT i 201311 then pre- 
dicts a black hole mass M = 1.6 x 10 5 M 0 . Given the 
intrinsic width of the M-L relation and the uncertainties 
in our black hole mass measurement, we conclude that, 
although a factor 3 higher, our mass measurement is not 
inconsistent with the M-L relation. 


6.4. Summary 

We presented dynamical models of the center of the 
nearest Seyfert 1 galaxy NGC 4395. NGC 4395 was the 
first bulgeless galaxy known to host an AGN and still 
has the most secure detection of an intermediate-mass 
black hole (M < 10 6 M 0 ), although with considerable 
discussion about its exact mass. 

Our high-resolution HST/WFC3 data enabled us to 
model the morphology and M/L ratio of the compact 
nuclear star cluster at the center of NGC 4395. After 
correcting for the asymmetric contribution of emission 
lines to the cluster light, we find that the cluster is al¬ 
most spherical (in projection). The range of M/L ratios 
of the cluster that we inferred from composite stellar pop¬ 
ulation models is in excellent agreement with those found 
for globular clusters in M31. 

Combining our HST imaging data with adaptive op¬ 
tics assisted near-IR integral field data of the cluster, 
we model the dynamics of the molecular gas. Our best 
fitting dynamical models contain an intermediate mass 
black hole with mass M = 4//® x 10 5 M 0 (3<r uncertain¬ 
ties). We characterize the influence of the uncertainty 
in the stellar potential on our measurement, and find 
that it leads to changes in the black hole mass by less 
than 1-er. Our mass estimate is in excellent agreement 
with th e mass estimate from reverberation mapping by 
iPeterson et al.l (|2005l) . but significantly highe r than the 
more recent estimate from lEdri et abl (2012) and esti¬ 
mates based on accretion. 

ACKNOWLEDGMENTS 

M.d.B. thanks Tom Maccarone for pointing out rel¬ 
evant references and Jenny Greene for comments on 
the manuscript. A.C.S. is supported by NSF CAREER 
grant AST-1350389. Research by A.,LB. is supported 
by NSF grant AST-1412693. M.C. acknowledges sup¬ 
port from a Royal Society University Research Fellow¬ 
ship. V.P.D is supported by STFC Consolidated grant 
# ST/J001341/1. Support for HST program GO-12163 
was provided by NASA through a grant from the Space 
Telescope Science Institute, which is operated by the 
Association of Universities for Research in Astronomy, 
Inc., under NASA contract NAS 5-26555. Based on ob¬ 
servations obtained at the Gemini Observatory, which 
is operated by the Association of Universities for Re¬ 
search in Astronomy, Inc., under a cooperative agree¬ 
ment with the NSF on behalf of the Gemini partnership: 
the National Science Foundation (United States), the Na¬ 
tional Research Council (Canada), CONICYT (Chile), 
the Australian Research Council (Australia), Ministerio 





























































The low mass MBH of NGC4395 


13 


da Ciencia, Tecnologia e Inovagao (Brazil) and Ministe- 
rio de Ciencia, Tecnologia e Innovation Productiva (Ar¬ 
gentina). Some of the data presented herein were ob¬ 
tained at the W.M. Keck Observatory, which is operated 
as a scientific partnership among the California Insti¬ 
tute of Technology, the University of California and the 
National Aeronautics and Space Administration. The 
Observatory was made possible by the generous finan¬ 
cial support of the W.M. Keck Foundation. The authors 
wish to recognize and acknowledge the very significant 
cultural role and reverence that the summit of Mauna 
Kea has always had within the indigenous Hawaiian com¬ 
munity. We are most fortunate to have the opportunity 
to conduct observations from this mountain. 

REFERENCES 


Abazajian, K. N., Adelman-McCarthy, J. K., Agiieros, M. A., 
Allam, S. S., Allende Prieto, C., An, D., Anderson, K. S. J., 
Anderson, S. F., Annis, J., Bahcall, N. A., Sz et al. 2009, ApJS, 
182, 543 

Antonini, F. 2013, ApJ, 763, 62 

Barth, A. J., Sarzi, M., Rix, H.-W., Ho, L. C., Filippenko, A. V., 
Sz Sargent, W. L. W. 2001, ApJ, 555, 685 
Barth, A. J., Strigari, L. E., Bentz, M. C., Greene, J. E., Sz Ho, 

L. C. 2009, ApJ, 690, 1031 

Begelman, M. C., Volonteri, M., Sz Rees, M. J. 2006, MNRAS, 
370, 289 

Bendinelli, O. 1991, ApJ, 366, 599 
Bertin, E. Sz Arnouts, S. 1996, A&AS, 117, 393 
Bertola, F., Cappellari, M., Funes, Jose G., S. J., Corsini, E. M., 
Pizzella, A., Sz Beltran, J. C. V. 1998, ApJ, 509, L93 
Boker, T., Schinnerer, E., Sz Lisenfeld, U. 2011, ASzA, 534, A12 
Bromm, V. 2013, Reports on Progress in Physics, 76, 112901 
Bromm, V. Sz Larson, R. B. 2004, ARA&A, 42, 79 
Cappellari, M. 2002, MNRAS, 333, 400 
—. 2008, MNRAS, 390, 71 

Cappellari, M., Verolme, E. K., van der Marel, R. P., Verdoes 
Kleijn, G. A., Illingworth, G. D., Franx, M., Carollo, C. M., Sz 
de Zeeuw, P. T. 2002, ApJ, 578, 787 
Carson, D. J., Barth, A. J., Seth, A. C., den Brok, M., 

Cappellari, M., Greene, J. E., Ho, L. C., Sz Neumayer, N. 2015, 
ArXiv e-prints 

Conroy, C., Gunn, J. E., Sz White, M. 2009, ApJ, 699, 486 
De Lorenzi, F., Hartmann, M., Debattista, V. P., Seth, A. C., Sz 
Gerhard, O. 2013, MNRAS, 429, 2974 
Desroches, L.-B., Filippenko, A. V., Kaspi, S., Laor, A., Maoz, 

D. , Ganeshalingam, M., Li, W., Moran, E. C., Swift, B., Bentz, 

M. C., Ho, L. C., Nandra, K., O’Neill, P. M., Sz Peterson, B. M. 
2006, ApJ, 650, 88 

Edri, H., Rafter, S. E., Chelouche, D., Kaspi, S., Sz Behar, E. 
2012, ApJ, 756, 73 

Emsellem, E., Monnet, G., Bacon, R., Sz, Nieto, J.-L. 1994, ASzA, 
285, 739 

Fan, X., Narayanan, V. K., Lupton, R. H., Strauss, M. A., Knapp, 
G. R., Becker, R. H., White, R. L., Pentericci, L., Leggett, 

S. K., Haiman, Z., Gunn, J. E., Ivezic, Z., Schneider, D. P., 
Anderson, S. F., Brinkmann, J., Bahcall, N. A., Connolly, 

A. J., Csabai, I., Doi, M., Fukugita, M., Geballe, T., Grebel, 

E. K., Harbeck, D., Hennessy, G., Lamb, D. Q., Miknaitis, G., 
Munn, J. A., Nichol, R., Okamura, S., Pier, J. R., Prada, F., 
Richards, G. T., Szalay, A., Sz York, D. G. 2001, AJ, 122, 2833 

Ferrarese, L. Sz Merritt, D. 2000, ApJ, 539, L9 
Filippenko, A. V. Sz Ho, L. C. 2003, ApJ, 588, L13 
Filippenko, A. V. Sz Sargent, W. L. W. 1989, ApJ, 342, Lll 
Fisher, D. B. Sz Drory, N. 2008, AJ, 136, 773 
Galleti, S., Federici, L., Bellazzini, M., Fusi Pecci, F., Sz Macrina, 
S. 2004, ASzA, 416, 917 

Gebhardt, K., Bender, R., Bower, G., Dressier, A., Faber, S. M., 
Filippenko, A. V., Green, R., Grillmair, C., Ho, L. C., 
Kormendy, J., Lauer, T. R., Magorrian, J., Pinkney, J., 
Richstone, D., Sz Tremaine, S. 2000, ApJ, 539, L13 


Gebhardt, K., Lauer, T. R., Kormendy, J., Pinkney, J., Bower, 

G. A., Green, R., Gull, T., Hutchings, J. B., Kaiser, M. E., 
Nelson, C. H., Richstone, D., Sz Weistrop, D. 2001, AJ, 122, 
2469 

Gebhardt, K., Richstone, D., Tremaine, S., Lauer, T. R., Bender, 
R., Bower, G., Dressier, A., Faber, S. M., Filippenko, A. V., 
Green, R., Grillmair, C., Ho, L. C., Kormendy, J., Magorrian, 
J., Sz Pinkney, J. 2003, ApJ, 583, 92 
Glover, S. 2005, Space Sci. Rev., 117, 445 
Graham, A. W. Sz Spitler, L. R. 2009, MNRAS, 397, 2148 
Greene, J. E. 2012, Nature Communications, 3 
Greene, J. E., Peng, C. Y., Kim, M., Kuo, C.-Y., Braatz, J. A., 
Impellizzeri, C. M. V., Condon, J. J., Lo, K. Y., Henkel, C., Sz 
Reid, M. J. 2010, ApJ, 721, 26 
Haehnelt, M. G. Sz Rees, M. J. 1993, MNRAS, 263, 168 
Hartmann, M., Debattista, V. P., Seth, A., Cappellari, M., Sz 
Quinn, T. R. 2011, MNRAS, 418, 2697 
Haynes, M. P., van Zee, L., Hogg, D. E., Roberts, M. S., Sz 
Maddalena, R. J. 1998, AJ, 115, 62 
Hicks, E. K. S., Davies, R. I., Maciejewski, W., Emsellem, E., 
Malkan, M. A., Dumas, G., Muller-Sanchez, F., Sz Rivers, A. 
2013, ApJ, 768, 107 

Hicks, E. K. S., Davies, R. I., Malkan, M. A., Genzel, R., Tacconi, 

L. J., Muller Sanchez, F., Sz Sternberg, A. 2009, ApJ, 696, 448 
Ho, L. C. Sz Kim, M. 2014, ApJ, 789, 17 

Jacobs, B. A., Rizzi, L., Tully, R. B., Shaya, E. J., Makarov, 

D. I., Sz Makarova, L. 2009, AJ, 138, 332 
Jedrzejewski, R. I. 1987, MNRAS, 226, 747 

Jiang, L., Fan, X., Bian, F., Annis, J., Chiu, K., Jester, S., Lin, 

H. , Lupton, R. H., Richards, G. T., Strauss, M. A., 
Malanushenko, V., Malanushenko, E., Sz Schneider, D. P. 2009, 
AJ, 138, 305 

Jordi, K., Grebel, E. K., Sz Ammon, K. 2006, A&A, 460, 339 
Karachentsev, I. D., Sharina, M. E., Dolphin, A. E., Grebel, 

E. K., Geisler, D., Guhathakurta, P., Hodge, P. W., 
Karachentseva, V. E., Sarajedini, A., Sz Seitzer, P. 2003, A&A, 
398, 467 

Kaspi, S., Smith, P. S., Netzer, H., Maoz, D., Jannuzi, B. T., Sz 
Giveon, U. 2000, ApJ, 533, 631 
Kormendy, J. Sz Ho, L. C. 2013, ARA&A, 51, 511 
Kormendy, J. Sz McClure, R. D. 1993, AJ, 105, 1793 
Kormendy, J. Sz Richstone, D. 1995, ARA&A, 33, 581 
Koushiappas, S. M., Bullock, J. S., Sz Dekel, A. 2004, MNRAS, 
354, 292 

Kraemer, S. B., Ho, L. C., Crenshaw, D. M., Shields, J. C., Sz 
Filippenko, A. V. 1999, ApJ, 520, 564 
Krajnovic, D., Cappellari, M., de Zeeuw, P. T., Sz Copin, Y. 2006, 
MNRAS, 366, 787 

Krist, J. E., Hook, R. N., Sz Stoehr, F. 2011, in Society of 
Photo-Optical Instrumentation Engineers (SPIE) Conference 
Series, Vol. 8127, Society of Photo-Optical Instrumentation 
Engineers (SPIE) Conference Series, 0 
Kroupa, P. 2001, MNRAS, 322, 231 
Kruijssen, J. M. D. Sz Mieske, S. 2009, A&A, 500, 785 
Lira, P., Lawrence, A., O’Brien, P., Johnson, R. A., Terlevich, R., 
Sz Bannister, N. 1999, MNRAS, 305, 109 
Lotz, J. M., Telford, R., Ferguson, H. C., Miller, B. W., Stiavelli, 

M. , Sz Mack, J. 2001, ApJ, 552, 572 

Macchetto, F., Marconi, A., Axon, D. J., Capetti, A., Sparks, W., 
Sz Crane, P. 1997, ApJ, 489, 579 
Magorrian, J., Tremaine, S., Richstone, D., Bender, R., Bower, 

G., Dressier, A., Faber, S. M., Gebhardt, K., Green, R., 
Grillmair, C., Kormendy, J., Sz Lauer, T. 1998, AJ, 115, 2285 
Matthews, L. D., van Driel, W., Sz Gallagher, III, J. S. 1998, AJ, 
116, 1169 

McConnell, N. J. Sz Ma, C.-P. 2013, ApJ, 764, 184 
McGregor, P. J., Hart, J., Conroy, P. G., Pfitzner, M. L., 
Bloxham, G. J., Jones, D. J., Downing, M. D., Dawson, M., 
Young, P., Jarnyk, M., Sz Van Harmelen, J. 2003, in Society of 
Photo-Optical Instrumentation Engineers (SPIE) Conference 
Series, Vol. 4841, Instrument Design and Performance for 
Optical/Infrared Ground-based Telescopes, ed. M. Iye Sz 
A. F. M. Moorwood, 1581-1591 

McHardy, I. M., Koerding, E., Knigge, C., Uttley, P., Sz Fender, 
R. P. 2006, Nature, 444, 730 



14 


den Brok, et al. 


McLean, I. S., Becklin, E. E., Bendiksen, O., Brims, G., Canfield, 
J., Figer, D. F., Graham, J. R., Hare, J., Lacayanga, F., Larkin, 
J. E., Larson, S. B., Levenson, N., Magnone, N., Teplitz, H., Sz 
Wong, W. 1998, in Society of Photo-Optical Instrumentation 
Engineers (SPIE) Conference Series, Vol. 3354, Infrared 
Astronomical Instrumentation, ed. A. M. Fowler, 566-578 
Merritt, D., Ferrarese, L., Sz Joseph, C. L. 2001, Science, 293, 

1116 

Mieske, S., Frank, M. J., Baumgardt, H., Liitzgendorf, N., 
Neumayer, N., Sc Hilker, M. 2013, A&A, 558, A14 
Minezaki, T., Yoshii, Y., Kobayashi, Y., Enya, K., Suganuma, M., 
Tomita, H., Koshida, S., Yamauchi, M., & Aoki, T. 2006, ApJ, 
643, L5 

Moran, E. C., Eracleous, M., Leighly, K. M., Chartas, G., 
Filippenko, A. V., Ho, L. C., Sz Blanco, P. R. 2005, AJ, 129, 
2108 

Moran, E. C., Filippenko, A. V., Ho, L. C., Shields, J. C., Belloni, 
T., Comastri, A., Snowden, S. L., Sz Sramek, R. A. 1999, 

PASP, 111, 801 

Nardini, E. Sz Risaliti, G. 2011, MNRAS, 417, 2571 
Neumayer, N., Cappellari, M., Reunanen, J., Rix, H.-W., van der 
Werf, P. P., de Zeeuw, P. T., Sc Davies, R. I. 2007, ApJ, 671, 
1329 

Neumayer, N. Sz Walcher, C. J. 2012, Advances in Astronomy, 
2012, 15 

Oiling, R. P. 1996, AJ, 112, 457 

Onken, C. A., Ferrarese, L., Merritt, D., Peterson, B. M., Pogge, 
R. W., Vestergaard, M., Sz Wandel, A. 2004, ApJ, 615, 645 
Peletier, R. F. Sc de Grijs, R. 1998, MNRAS, 300, L3 
Peng, C. Y., Ho, L. C., Impey, C. D., Sc Rix, H.-W. 2002, AJ, 

124, 266 

Peterson, B. M., Bentz, M. C., Desroches, L.-B., Filippenko, 

A. V., Ho, L. C., Kaspi, S., Laor, A., Maoz, D., Moran, E. C., 
Pogge, R. W., Sc Quillen, A. C. 2005, ApJ, 632, 799 
Riffel, R. A., Storchi-Bergmann, T., Dors, O. L., Sz Winge, C. 
2009, MNRAS, 393, 783 

Rossa, J., van der Marel, R. P., Boker, T., Gerssen, J., Ho, L. C., 
Rix, H.-W., Shields, J. C., Sc Walcher, C.-J. 2006, AJ, 132, 1074 
Sandage, A. Sz Tammann, G. A. 1981, A revised Shapley-Ames 
Catalog of bright galaxies 

Sarzi, M., Rix, H.-W., Shields, J. C., Rudnick, G., Ho, L. C., 
McIntosh, D. H., Filippenko, A. V., Sz Sargent, W. L. W. 2001, 
ApJ, 550, 65 

Schlafly, E. F. Sc Finkbeiner, D. P. 2011, ApJ, 737, 103 


Sersic, J. L. 1968, Atlas de galaxias australes 
Seth, A., Agiieros, M., Lee, D., Sc Basu-Zych, A. 2008, ApJ, 678, 
116 

Seth, A. C., Cappellari, M., Neumayer, N., Caldwell, N., Bastian, 
N., Olsen, K., Blum, R. D., Debattista, V. P., McDermid, R., 
Puzia, T., Sc Stephens, A. 2010, ApJ, 714, 713 
Seth, A. C., Dalcanton, J. J., Hodge, P. W., Sz Debattista, V. P. 
2006, AJ, 132, 2539 

Seth, A. C., van den Bosch, R., Mieske, S., Baumgardt, H., den 
Brok, M., Strader, J., Neumayer, N., Chilingarian, I., Hilker, 
M., McDermid, R., Spitler, L., Brodie, J., Frank, M. J., Sz 
Walsh, J. L. 2014, Nature, 513, 398 
Shapiro, K. L., Cappellari, M., de Zeeuw, T., McDermid, R. M., 
Gebhardt, K., van den Bosch, R. C. E., Sz Statler, T. S. 2006, 
MNRAS, 370, 559 

Spinelli, P. F., Storchi-Bergmann, T., Brandt, C. H., Sz Calzetti, 
D. 2006, ApJS, 166, 498 

Strader, J., Caldwell, N., Sz Seth, A. C. 2011, AJ, 142, 8 
Strader, J., Smith, G. H., Larsen, S., Brodie, J. P., Sz Huchra, 

J. P. 2009, AJ, 138, 547 

Swaters, R. A., Schoenmakers, R. H. M., Sancisi, R., Sc van 
Albada, T. S. 1999, MNRAS, 304, 330 
Thim, F., Hoessel, J. G., Saha, A., Claver, J., Dolphin, A., Sz 
Tammann, G. A. 2004, AJ, 127, 2322 
van der Marel, R. P. Sc van den Bosch, F. C. 1998, AJ, 116, 2220 
Vaughan, S., Iwasawa, K., Fabian, A. C., Sc Hayashida, K. 2005, 
MNRAS, 356, 524 

Venemans, B. P., Findlay, J. R., Sutherland, W. J., De Rosa, G., 
McMahon, R. G., Simcoe, R., Gonzalez-Solares, E. A., Kuijken, 

K. , Sc Lewis, J. R. 2013, ApJ, 779, 24 

Verdoes Kleijn, G. A., van der Marel, R. P., Carollo, C. M., Sz de 
Zeeuw, P. T. 2000, AJ, 120, 1221 
Verolme, E. K., Cappellari, M., Copin, Y., van der Marel, R. P., 
Bacon, R., Bureau, M., Davies, R. L., Miller, B. M., Sc de 
Zeeuw, P. T. 2002, MNRAS, 335, 517 
Volonteri, M., Lodato, G., Sc Natarajan, P. 2008, MNRAS, 383, 
1079 

Walcher, C. J., Boker, T., Chariot, S., Ho, L. C., Rix, H.-W., 
Rossa, J., Shields, J. C., Sz van der Marel, R. P. 2006, ApJ, 

649, 692 

Willott, C. J., Delorme, P., Reyle, C., Albert, L., Bergeron, J., 
Crampton, D., Delfosse, X., Forveille, T., Hutchings, J. B., 
McLure, R. J., Omont, A., Sz Schade, D. 2009, AJ, 137, 3541 
Wrobel, J. M. Sz Ho, L. C. 2006, ApJ, 646, L95 



The low mass MBH of NGC4395 


15 



2.290 2.295 2.300 2.305 2.310 

Rest wavelength [gm] 


Figure 9. Normalized NIRSPEC spectra of the NGC 4395 nucleus and a K3III star. The wavelength of the first CO bandhead is 2.2935 
/zm. We also show the expected spectrum of the nucleus based on modeling of the K3III star for 65% (blue dashed line) and 95% AGN 
contribution (red dotted line) and a velocity dispersion of a = 30 kms -1 . The AGN is clearly dom inating the spectrum. This domination 
of the continuum component is also easy to discern from the Gemini/NIFS data shown in Fig. llOl where even in the least AGN dominated 
aperture, the line depths of the CO bandheads are at most a few percent. 


APPENDIX 

A. STELLAR VELOCITY DISPERSION: NIRSPEC & NIFS DATA 

The stellar velocity dispersion of the nuclear cluster in NGC 4395 is a para meter of interest for dynamica l modeling 
of th e cluster providing a constraint on the total mass of the system (e.g. iBarth et al.l 120091 iNeumaver fc Walcheil 

EoU). 

To date, the on ly direct constraint on the velocity dispersion of this cluster is a measurement of <r < 30 kms” 1 by 
iFilippenko fc Hoi (I2003T) based on high-dispersion optical spectroscopy around the Ca II triplet obtained with HIRES 
at the Keck Observator y. As an attempt t o measure a directly, we obtained near-infrared observations with the 
NIRSPEC spectrograph (IMcLean et al.lfl998l) at the Keck II telescope. 

NIRSPEC observations were obtained during the first half of two nights, 2010 March 28-29 UT, using the NIRSPEC 
echelle mode with a 0"72 x 12" slit. Weather was clear and the seeing was excellent, with focus runs giving FWHM 
of ~ 0"4 on stellar images. On the first night we observed NGC 4395 in the NIRSPEC-7 (AT-band) filter, and during 
the second night the NIRSPEC-5 (F/-band) filter was used. The galaxy nucleus was nodded along the slit between 
successive exposures for sky subtraction, with individual exposure times of 300 or 600 s, and the slit was held at a 
fixed position angle of 120°. The total on-source exposure time was 2.0 hr in the K band and 2.3 hr in the H band. 
We observed A0V stars as telluric calibrators, and several K and M giant stars as velocity templates. Arc lamp and 
flat-field exposures were obtained at the start and end of each half-night. 

Reductions were carried out using the NIRSPEC REDSPEC pipeline. For the A'-band data, we focused on the 
spectral order containing the primary CO bandhead at 2.2935 /im since this is the main feature of interest for measuring 
a. In this spectral order, the dispersion is 3.33 x 10” 5 n m pixel” 1 . Extracted spectra were wavelength calibrated and 
flux calibrated using arc lamp and A0V star exposures, and individual exposures were combined to give a weighted 
average spectrum. Figure [9] shows the CO bandhead region spectrum of the NGC 4395 nuclear cluster and of a K3III 
giant for comparison. Over this region, the NGC 4395 spectrum has S/N = 23 pixel” 1 . There is no detected CO 
bandhead feature in NGC 4395. If 95% of the spectrum were coming from the AGN, we would not be able to detect 
CO features at a velocity dispersion of 30 kms” 1 . If the AGN contributes only 65% of the flux, the features would 
be clearly visible. Apparently the near-IR spectrum of the NGC 4395 nucleus is so dominated by the AGN (and 
possibly by a contribution from young, hot stars) that the CO features are diluted to the point of being invisible in the 
integrated light of the cluster and AGN. This domination of the AGN continuum in the near IR is supported by the 
ground-based near-IR observations of iMinezaki et all (12006!) . which show clear variability in NGC 4395’s integrated 
nuclear flux by a factor > 2. 

The H -band spectrum is significantly noisier and also does not show any stellar absorption features from the cluster. 
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Figure 10. Normalized NIFS spectrum of the NGC 4395 nucleus. The first two CO bandheads are marked by the dashed lines. Also 
marked (by dotted lines) are the H 2 3-2 S( 2) line at 2.28703 //rri and the Ca VIII line at 2.32204 //rn. 


We also searched our NIFS data for signs of the CO bandhead. We tested various annuli for signs of the CO 
bandheads. Annular spectra with radii between 0"15 and 0"35 show the strongest signs of the first CO bandhead. 
Only the first CO bandhead is seen at a depth of 3%, a nd the line is at a redshift of 317 kins -1 , within a couple 
kms -1 of the redshift of the galaxy derived from H I gas (lHavnes et al.lfl99 8f). The formal fit to the dispersion of the 
system suggests a dispersion of the system <30 kms -1 ; however due the weakness of the line, the lack of additional 
CO bandhead features and Ca VIII emission line contamination of the 2nd CO ban dhe ad (see Fig. HOl) . we regard this 
measurement as less reliable than the Calcium triplet measurement shown in lFilippenko fe Hoi (120031) . 

It would still be extremely desirable to measure the velocity dispersion of this cluster, but it is probably necessary 
to couple such a measurement with monitoring of the flux level of the AGN to ensure minimal dilution of the spectral 
features by the AGN component. 
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Figure 11. Results of dynamical models while using photometric model B (iterative residual subtraction). The plots show the likelihood 
of the data for models with free inclination, M/L ratio and black hole mass plotted Left: as a function of M/L ratio and black hole mass, 
maximized for inclination. Right : as a function of inclination and black hole mass, maximized for M/L. Contours encompass models 
within 1 and 3 cr distance from the best-fit model. The best fit model is highly inclined. 
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Figure 12. Same as Fig. 1111 but now for the star cluster model derived after subtracting the F275W image. A high black hole mass is 
preferred. 


B. DYNAMICAL MODELING RESULTS FOR ALTERNATIVE MODELS OF THE CLUSTER 

A major uncertainty in our modeling of the gas kinematics is the exact structure and mass of the stellar potential, 
which is difficult to measure both because of the NSC’s compactness and because of the presence of line emission. In 
Section O we presented, besides ordinary masking of the emission line regions, two alternative ways to deal with them. 
The different estimates for the M/L ratios of the stellar populations for these alternative models are presented in Fig. 
[4j Here we present also the results of the dynamical models for these alternative models. Fig. [TT| shows the results 
of the dynamical M/L and inclination versus the MBH mass for the second structural parameter model (iterative 
residual subtraction). Despite the 10% difference in NSC mass and different Sersic index of the profile, the black hole 
mass is in agreement with the one from the analysis in the main text. The modeling results using the other structural 
parameter model are shown in Fig. [T2] The best-fit BH mass for this model is higher than the other two models, and 
si mple axisymmetri c Jean s models show that there is some tension with the stellar velocity dispersion measurement 
of lFilinnenko fe Hoi (12003!) . Possible explanations are that the F275W band contained not just line emission but also 
stellar emission, or that the scaling between the F275W band and the other bands is not the same for NGC 4395 and 
NGC 1068. 
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Figure 13. Photometric decomposition of the SDSS i band data of NGC 4395 into an exponential and Sersic model. The central Sersic 
model has an effective radius of 26 // (0.5 kpc), and a Sersic index n = 0.43. 

C. PHOTOMETRIC EVIDENCE FOR A BULGE IN NGC 4395 

As far as we are aware, there exists no photometric decomposition of NGC 4395 into a bulge and a disk. Such a 
decomposition is important to corroborate the claim that MBHs exist in bulgeless galaxies. We therefore present a 
photometric decomposition of this galaxy here. 

We use i band imaging data from the Sloan Digital Sky Survey (SDSS, lAbazaiian et al.l 12009 1 to constrain the 
presence or absence of a bulge in NGC 4395. The data from this band are the least affected by dust and young stellar 
populations and less shallow than the 2 -band data. We note that NGC 4395 falls close to the edge of the chip. This 
does not affect our analysis in the inner parts, but beyond 220" the data are dominated by a single irregular-looking 
spiral arm because of thi s asymmetry. __ 

We use SExtractor (|Bertin fe ArnoutsH TOOO ) to identify background galaxies, for eground stars and star clusters 
in the image and generate a bad pixel mask. We then use the IRAI0 task ellipse (|.Tedrzeie wski 1983) to measure 
the radial surface brightness profile of the galaxy. We fit the profile with the combination of an exponential profile 
and a Sersic profile. We exclude the nucleus from the fit, as well as the data beyond 220". 

Fig. US shows the decomposition of NGC 4395 into a Sersic and an exponential model. The exponential model has 
magnitude rip = 10.3 mag (Mj = —17.9 mag), scale length h = 102" (2.2 kpc) and dominates over the Sersic profile at 
all radii. The Sersic component has up = 13.7 mag (Mj = —14.5 mag), effective radius R e s = 26" (0.5 kpc) and Sersic 
index n = 0.44. The low Sersic i ndex and low luminosity (the bulge-to-total ratio B/T = 0.04) indicate that this is 
probably not a classical bulge IlFisher fc Drorvll2008 ). but more likely a bar. Visually, the elongated shape suggests 
that it is indeed a bar, and we also find tentative support for this by a small change of 10 degrees in the position angle 
in the inner 40 arcsec of NGC 4395. 


11 IRAF is distributed by the National Optical Astronomy Ob¬ 
servatory, which is operated by the Association of Universities for 


Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 













